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We show that all so{N)i universality class quantum criticalities emerge when one-dimensional gen¬ 
eralized cluster models are perturbed with Ising or Zeeman terms. Each critical point is described by 
a low-energy theory of N linearly dispersing fermions, whose spectrum we show to precisely match 
the prediction by so{N)i conformal field theory. Furthermore, by an explicit construction we show 
that all the cluster models are dual to non-locally coupled transverse field Ising chains, with the 
universality of the so(A)i criticality manifesting itself as N of these chains becoming critical. This 
duality also reveals that the symmetry protection of cluster models arises from the underlying Ising 
symmetries and it enables the identification of local representations for the primary fields of the 
so{N)i conformal field theories. For the simplest and experimentally most realistic case that corre¬ 
sponds to the original one-dimensional cluster model with local three-spin interactions, our results 
show that the su{2)2 — so(3)i Wess-Zumino-Witten model can emerge in a local, translationally 
invariant and Jordan-Wigner solvable spin-1/2 model. 


A striking property of quantum phase transitions is 
that of universality. Near a critical regime separating 
distinct quantum states of matter, the system specific 
microscopic details are lost and the system acquires uni¬ 
versal behavior that is characterized by scaling exponents 
of distinct observables near the transition [1]. A crucial 
step in putting order to the zoo of distinct universali¬ 
ties was made by noting that the diverging correlation 
length at the critical point implies conformal invariance, 
and thus a description by a conformal field theory (CFT) 
[HE]. This approach is particularly powerful in one spa¬ 
tial dimension (ID), where the relevant CFT fully de¬ 
scribes (the scaling of) the correlation functions. 

ID quantum models with critical points are thus the 
natural playgrounds for quantum criticality. The sim¬ 
plest and the most celebrated ones are the transverse field 
Ising (TFI) chain [3] , that at criticality is described by the 
Ising CFT, and the XY chain where the anisotropic 
transition is described by the so called so(2)i CFT. So 
ubiquitous are these universality classes that criticality 
beyond them is often filed under “exotic”. Still, consid¬ 
ering this ubiquitous nature these simple models have 
played in various fields of physics, it was only rather 
recently when the details of the criticality of the TFI 
chain were experimentally probed [B]. The challenge of 
probing quantum criticality beyond these simple models 
owes to the lack of tractable models that admit acces¬ 
sible experimental realization. To go beyond them, of¬ 
ten higher spins [THU], infinite range couplings Haul, 
broken translational invariance HU, strongly interacting 
systems [TMTH] - or combinations there of [B HD] , are 
required. Here we show that a hierarchy of local, ex¬ 
actly solvable and translationally invariant spin-1 /2 mod¬ 
els provides the simplest setting to generalize the Ising 


universality class by realizing all so{N)i criticalities, that 
were previously considered in manifestly so{N) symmet¬ 
ric higher spin chains [3 Ej or in models with broken 
translational invariance m- Furthermore, we show that 
the universality manifests itself microscopically through 
every quantum critical point being dual to N non-locally 
coupled critical transverse field Ising chains. 

Our hierarchy can be viewed as a generalization of ID 
cluster models [H] - a class of stabilizer models with 
a symmetry protected degenerate ground state manifold 
[B ED] - whose ground states in 2D were first proposed 
as a universal resource for one-way quantum computa¬ 
tion El]. While in ID these models are universal re¬ 
sources only for a single qubit, they have proven acces¬ 
sible settings to study the entanglement E5H5D] 

and the computational power I31HTI) of symmetry pro¬ 
tected states, as well as the robustness of edge states in 
a many-body localized phase [34]. Perturbing the pure 
stabilizer models with Ising and Zeeman terms, we define 
the generalized cluster models with periodic boundary 
conditions by the Hamiltonians (the original ID cluster 
model corresponds to A = 3) 

L-l 

^iluLer = E ( + •^«+l + ) > (1) 

where Cf = crfcrf+i • • • <+(^_ 2 ) 0 -f+(^_i) acts on A adja¬ 
cent sites and J and h are the magnitudes of the Ising 
and Zeeman terms, respectively. The cluster stabilisers 
Cf commute with each other and thus for J = h = 0 
one obtains the cluster state as the unique ground state 
that satisfies Cf | 4') = — | 4/) for all i. Anologous to 
the spin-1 Haldane phase [DDHDZ], these states are sym¬ 
metry protected topological (SPT) phases protected by 
global symmetries = Hiflo c^fA-l)^+n 
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FIG. 1. Phase diagram of the 3-cluster model exhibiting 
gapped cluster (C), polarised (P) and anti-ferromagnetic 
{AFM^''^) phases |27| . The solid (dashed) lines indicate sec¬ 
ond order quantum phase transitions along which the spec¬ 
tral gap closes at two (one) distinct Fermi momenta. The 
insets show the spectrum ® at the multi-critical points 
(J, h) = (±1,0) and (J, h) = (0,±1) where the critical spec¬ 
trum exhibits two (circles) and three (squares) Majorana 
cones, respectively. 


for n = 0,...,A—2, and characterized by a string order 
parameter and edge states HSl |3H] ■ We show below 
that even if these symmetries are broken, the perturbed 
Hamiltonians can still exhibit global symmetries, 
where N depends on the perturbation in question. 

To study the phase diagrams of the H-cluster models 


(121, we map (121 to a problem of free fermions by intro¬ 


ducing two Majorana fermions per physical site 



h=(l[> 


( 2 ) 


\k<i 


Assuming a chain of L sites with periodic boundary 
conditions, in momentum space the Hamiltonian takes 


where Ck = 


the form = iJ2k 

—h + Jexp{—ik) + exp[j(A — l)fc] and the Majorana op¬ 
erators satisfy in the momentum space c\. = C-k- The 
momentum takes the values k = (2m ± 1 — P)'k/L, 
where m = 0, ...,L — 1 and P = ±1 is the eigen¬ 
value of the parity operator Vduster = Il^o^ 
spectrum can be brought to the diagonal form H = 
Sfc kfc|( 27 fe 7 fc ~ 1) by moving to the complex fermion 

basis Ik = :^ + bk^ with the eigenvalues given 

by 

|efc| = [l ± ± ± 2 J cos(Afc) 

—2h{J cos(/c) + cos((A — l)fc))]^^^. (3) 


The resulting phase diagram for the simplest 3-cluster 
model, that was first studied in Ref is shown in 
Figure While the phase diagrams become in general 
more complex for A > 3 (presented in the Supplementary 


Material [55]), they all share a common feature: One al¬ 
ways finds multi-critical points at {J,h) = (±1,0) and 
(J, h) = (0, ±1), where there is a second order quan¬ 
tum phase transition from the cluster phase to an anti¬ 
ferromagnetic or a spin polarised phase, respectively. De¬ 
termining the CFT describing these critical points for 
general A is our main result. 

To gain insight into these critical points, we focus first 
on (J, h) = (±1, 0) where the dispersion reduces to |e/j| = 
v/2^ cos(Afc). This vanishes at the N = A Fermi 
points Kn = n/N + 2Trn/N (J = 1) or Kn = 2Tm/N 
(J = —1), where n = 0,..., iV — 1. Expanding around 
each Fermi momenta by writing k = Kn +p with p 1, 
one finds linear dispersion |eic^+p| = Np + 0{p^) im¬ 
plying a low-energy description in terms of N fermions 
of velocity N. Carrying out a similar analysis for the 
(J, h) = (0,±1) critical points, one finds N = A — 1 
Fermi points with linearly vanishing dispersion, as illus¬ 
trated in Figure [l] It is well known that at low energy, 
the spectrum of the critical TFI chain is described by a 
single linearly dispersing fermion and that the critical¬ 
ity is described by the Ising CFT with central charge 
c = 1/2 12 . Thus the low-energy description in terms 
of N fermions with a linear dispersion naively suggests 
that the corresponding critical points are described by 
a product theory of N Ising CFTs with c = N/2. In¬ 
deed, central charges of c = 1 and c = 3/2 have been 
obtained for the 3-cluster model at {J,h) = (0,±1) [5U] 
and {J,h) = (±1,0) |23|, respectively. 

However, neither the number of linearly dispersing 
fermions nor the central charge uniquely fix the CFT. In 
addition to the product CFT Ising^'^, the central charge 
of c = N/2 and the spectrum of N fermions are also 
consistent with the so called so(A^)i CFTs [7II51 [TK]. As 
these theories have dramatically different primary field 
content |51|55|, they can be distinguished by the energy 
levels and degeneracies in the finite-size energy spectrum 
that is fully determined by the field content of the corre¬ 
sponding CFT. For non-chiral models such as ours, the 
spectrum of an L site chain takes the form |3| 


E — EqL — 


TTVC 



2'kv 


{2ha + n), 


(4) 


where the on-site energy Eq and the velocity v are non- 
universal numbers and n is a non-negative integer. On 
the other hand, the central charge c and the scaling di¬ 
mension ha of each primary field a are universal and 
determined by the CFT. The spectrum can be simplified 
by setting the ground state energy to A = 0 and scaling 
the spectrum such that the first excited state has energy 
E = 2ha '1 with ha' being the smallest non-zero scal¬ 
ing dimension. The spectrum assumes then the simple 
form E = 2ha+n where the energies arrange themselves 
into integers offset by the scaling dimensions that char¬ 
acterise the CFT. The predictions of the energy levels 
and degeneracies by both so{N)i and Ising^-^ CFTs are 
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A=3 cluster model A=3 cluster model 

L=:36, J=l, h=0, so(3)^ point L=32, J=0, h=-l, so(2)j point 



where Qf = U^<jiTiiTii+ 2 ) and Q> = ]!*>:, 

Assuming that L is a multiple of three, applied to (12) 
we obtain 


H. 


(3) 


cluster 


L/3-1 

= 


X X 

3j^3j+3 


+T3^ + JViV2rl_^rS + 


j=0 


+Jri^+i + VoV2Tl_2rf+ 
T^j+2T^,+5 + Jt|,+2 + VoVitI_^t^) , 

( 6 ) 


where the duals of Vn = Y\m=o^ '''sm+n are hidden global 
symmetries of the perturbed cluster model (all the 
Zn symmetries are broken by J ^ 0). A simpler duality 
exists for h ^ 0, but J = 0 when the symmetries are 
preserved. Assuming that L is even, we introduce 


_ n<r*' 

^2j - ^2j> 


_ n>ry 

^2j + l - ^2i+H 


(7) 


FIG. 2. Matching of the CFT spectra E = 2ha + n for the 
so(3)i and so(2)i critical points of the N = 3 cluster model. 
The so(3)i CFT has three primary fields {1, ip, a} with scaling 
dimensions hi — 0, h^j, = 1/2 and ha = 3/16, while the so(2)i 
CFT has four primary fields A,A} with scaling dimen¬ 

sions hi = 0, h,p = 1/2 and hx = h^ = 1/8. In both cases the 
rescaled energies follow the predicted pattern Q, with the 
degeneracies (not indicated) at each energy level matching 
those predicted for so(A)i CFTs, see |15| . 


presented in the Supplementary Material [39]. In Fig¬ 
ure we plot the spectra for the 3-cluster model at the 
critical points (J, h) = (1,0) and (0,1), which match pre¬ 
cisely the prediction by the so(3)i and so(2)i CFTs, re¬ 
spectively. We have verified that this structure holds for 
general A-cluster models, which leads to our main result: 
The multi-critical points {J,h) = (±1,0) and (±1,0) of 
every A-cluster model have low-energy theories in terms 
N = A and N = A — 1 fermions, respectively, and they 
are described by so(A^)i CFTs. 

To analytically demonstrate the universality of these 
critical points, we now turn to show that at each so{N)i 
critical point the system can be mapped to N non-locally 
coupled TFI chains with symmetry. These symme¬ 
tries correspond to different decompositions of the parity 
sectors Vduster = Yltzl 2.n = 11^=1 'Pn i'^to hidden sym¬ 
metry sectors labeled by Vn = ±1- The canonical spin 
duality transformations we employ are inspired by Refs 
[13 im and exist for general N (presented in the Supple¬ 
mentary Material [33). For the sake of clarity we focus 
here on the 3-cluster model. When h = 0 we employ first 
the canonical transformations 


cr 

a 


4 

V 

3j + l 

V 

33+2 


— 0<T^ 

- ±^ + 1) 


'3i-l-2 


Q/ 


_ V V Q> 

^33 - ^3V 33 + 1^3 ’ 

Z _ ^<—Z X X 

^33 + 1 - ^3 ^33^3j + lT^3j+2 
^33+2 — T^3V3j+2^ 


(5) 

~^3 


where now Qf = t^i+i and Q> = lli>j ^ 21 , which 
directly give 

L/2-1 

tt{ 3) _ +hT^ T^-l- 

cluster ~ \j 2 V 23+2 + ' l^L-2^0 + 

3=0 

T-fj + lTfj+3 +hT^j + i+roTl_,^Tf, 

( 8 ) 

where Vn — '^ 2 m+n aow coincide with the Z„ 

symmetries. 

Along both cuts of the phase diagram that contain the 
multi-critical points we thus find the same structure: The 
A-cluster models with Ising or Zeeman perturbations are 
dual to TFI chains that are non-locally coupled such that 
the boundary condition for the nth TFI chain depends 
on the product Y[m=^n Pm of fh® Ising symmetry sectors 
of all the other chains. This coupling is precisely of the 
general form that, based on the framework of condensate- 
induced transitions between topological phases [ID], has 
been recently shown to change the criticality between 
Ising^^ and so(A^)i universality classes [T3]. In the sec¬ 
tor where Vn = I for all n, or for open chains with free 
boundary conditions, the spectrum is indeed equal to de¬ 
coupled TFI chains and at criticality the spectrum is 
indistinguishable from that of Ising^-^ CFT. The cru¬ 
cial point is that this sector is only a single subsector 
of the Vduster = 1 symmetry sector. To construct the 
full so{N)i critical spectrum, one needs to appropriately 
combine the spectra of N critical TFI chains using all 
Ising symmetry sectors. This subtle mixing of the spectra 
from different symmetry sectors underlies the emergence 
of the so{N)i instead of the Ising universality class 
and highlights the need to consider all symmetry sectors 
when identifying the relevant CFT by spectral means. 

These results can be generalized to the critical points 
between arbitrary cluster phases. Consider the critical 
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(J.h)=(±1.0) 



FIG. 3. The Majorana representation of the clus¬ 

ter model. The ellipses denote the physical spin sites that 
are represented by two Majoranas. The solid lines corre¬ 
spond to Majorana terms originating from afaf^i = ibiai+i 
or af = —iaibi terms, while dashed lines correspond to clus¬ 
ter terms = iaibi+ 2 - Along the h = 0 (J = 0) 

cut the system decouples locally into three (two) Majorana 
wires. The wires are coupled non-locally though through their 
boundary conditions that for each chain are either periodic or 
anti-periodic depending on Vduster = Uji'ia.jbj) = ±1. 

Hamiltonian H = where we have in¬ 

troduced two non-commuting cluster stabilizers = 

’^ith a = x,y (we define 
j = —erf ). Employing the mapping to Majorana 
fermions (i), one immediately finds that for a ^ j3 the 
system maps into |(H — 1) -|- (H — 1)| Majorana chains, 
which are only coupled through their boundary condi¬ 
tions given by Vciuster (see Figurej^for the Majorana rep¬ 
resentation of the 3-cluster model with A = 3 and B = 1 
or 2). For a = f3 one hnds |H —i?| such chains. Since Ma¬ 
jorana chains are well known to be dual to TFI chains, 
their coupling through their boundary conditions trans¬ 
lates precisely into the non-local coupling between TFI 
chains discussed above. Thus we conclude that the criti¬ 
cal points between arbitrary cluster phases are always de¬ 
scribed by the so{N) i CFT with N = |(yl—l)±(i?—1)|. 
That all critical points due to arbitrary competing cluster 
stabilizers can be mapped to the Hamiltonian 

L-N N L 

HsoiN). = E 

i—l n—1 m^n i—1 

(9) 

with Z 2 ^ symmetry, is the microscopic manifestation of 
the universality of so(7V)i criticality and provides a full 
classiheation of criticalities in this set of exactly solvable 
models. 

This universal microscopic description can be em¬ 
ployed to identify candidates for the spin representa¬ 
tions for fields related to the chiral primaries if and a 
of the so(A^)i CFT with scaling dimensions = 1/2 
and ha = iV/16, respectively. If a field with scaling di¬ 
mension Aq, = 2ha is represented by a local operator Of , 
then the ground state correlator at the critical point is 
expected to decay as {OfOf) ~ For a criti¬ 

cal TFI chain, which can be viewed as the special N = \ 


case of the hierarchy, the non-chiral combinations of the 
Ip = 1/2) and a {ha = 1/16) primaries are related 
to local operators through ipLipR ^ and aL<JR ^ 
(the chiral left (L) and right (R) moving parts require 
non-local operators) [5] . Thus local operators in the A- 
cluster models with the scaling h = 1 of the ipL'pR field 
are given by the duals of the rf operators, while a natu¬ 
ral candidate for a local operator with scaling dimension 
27V/16 of (JR(JR is the dual of the product of N adjacent 
Tj operators. For the so(3)i critical point of the 3-cluster 
model at (J, h) = (±1,0) these would explicitly be given 
by, for instance, ipL'pR ~ and (JL<yR cA. We 

leave the explicit verification, either analytically [4l [5] or 
numerically [301 mmi for future work. 

Finally, we show that the dual picture in terms of N 
boundary coupled TFI chains is consistent with and use¬ 
ful also to study the perturbed gapped cluster phases 
mM\- For J = h = 0 and for periodic boundary condi¬ 
tions, for A odd (even) there holds 11^1 = IliLi 

(1), which implies L (L—l) independent constraints and 
thus a unique (two-fold degenerate) ground state. In the 
cluster phase A — 1 TFI chains are in the ferromagnetic 
phase, but the unique ground state for A odd arises from 
the boundary coupling enforcing that there is only one 
sector {Vn = 1 for all n) where all chains have periodic 
boundary conditions (the lowest energy contribution per 
TFI chain), while for A even this occurs for two sec¬ 
tors {Vn = ±1 for all n). For free boundary conditions 
the boundary coupling is removed and there are only 
L—{A—1) stabiliser constraints. This gives rise to a 
fold ground state degeneracy that is consistent with A—1 
out of the N unconstrained TFI chains being in the fer¬ 
romagnetic phase. While the Zn symmetries protecting 
the pure cluster state are in general broken by Zeeman, 
Ising or competing cluster term perturbations, the dual 
TFI chain picture with hidden symmetry persists 

throughout the cluster phase. Since this picture enables 
to isolate the hidden degrees of freedom that contribute 
to the properties of perturbed cluster states (regardless 
of N only A—1 TFI chains contribute to the H-cluster 
state), it can be valuable for studying the computational 
power of symmetry protected topological states |5TH55] . 

We have shown that transitions between symmetry 
protected cluster phases realize all quantum critical 
points in the universality class of so{N)i CFTs. We 
explicitly demonstrated the microscopic universality of 
the transitions by mapping every multi-critical point to 
N non-locally coupled TFI chains with Zf^ symmetry. 
Thus in addition to being of interest to quantum infor¬ 
mation , cluster models are also accessible 

platforms for probing quantum criticality beyond Ising 
universality class that can emerge between distinct SPT 
states. Experimental realizations of these models has 
been proposed in optical lattices [25] and with trapped 
ions |43] . Furthermore, an important corollary of our 
results is the discovery that the su{2)2 — so(3)i critical- 
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ity, which has previously been discovered only in spin-1 
or higher systems III [IS EH [20] - can emerge also in a 
translationally invariant spin-1/2 chain with local inter¬ 
actions. A full classification of criticalities between SPT 
phases beyond the protection of the integrable cluster 
models remains an interesting open question. 
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Supplementary Material 

In the supplementary material we present the phase 
diagrams and give the duality transformations for general 
^-cluster models. We also briefly review how the CFT 
predicts the spectrum at a critical point and give the 
predictions by so(iV)i and Ising^^ CFTs. 


Phase diagrams for arbitrary ^-cluster models 


The spectrum of the A-cluster model is in general given 

by 

|efc| = [l + + 2 J cos(Afc) 

—2h{Jcos{k) + cos((A — . (10) 

To find the parameters J and h for which the gap closes, 
we complete squares to obtain \ek\^ = {h — Jcos(fc) — 
cos((A — l)fc))^ + (jsin(A:) — sin((A — l)fc))^. We find 
that there is gap closing point aX k = 0 when J = h — 1 
and at fc = TT when J = In addition, the gap 

closes along a curve in (h, J), that can be parametrized 
by the momentum 0 < p < tt for which Cp = 0. This 
curve is given by 


{h,J){p) 


sin(>lp) 

sin(p) 


sin((A — l)p) 
sin(p) 


)■ 


( 11 ) 


Along this curve, the gap closes at both k = p and 
k = 2Tr — p. The multi-critical points at (h, J) = (±1, 0) 
and (h, J) = (0,±1) occur when this line intersects it¬ 
self and/or the two lines with a gap closing at a hxed 
momentum. 

The phase diagrams for some of the smallest A are 
illustrated in Figure As A increases, the A-cluster 
phase becomes sharply bounded by J -I- h < 1. Several 
fine-tuned gapped phases also emerge due to the compe¬ 
tition between all the three Hamiltonian terms. These 
are magnetic phases similar to the incommensurate anti¬ 
ferromagnetic phases of the 3-cluster model studied in 
Ref [27]. 


Duality transformations for arbitrary A-cluster 
models 


so{N)i universality class. With the aid of duality trans¬ 
formations a hierarchy of local, but staggered spin models 
with these criticalities were constructed. Here we show 
that this scheme can be simplified considerably by using 
another set of duality transformations that gives directly 
the translationally invariant Wcluster models considered 
in the main text. 

The appropriately boundary coupled system of N crit¬ 
ical TFI chains is described by the Hamiltonian 


H. 


(AT) 

cTFI 


N-1 

= E 

n=0 


LIN-2 


y ^ {jN j+n'^N{j+l)+n A 


j=0 




( 12 ) 


where Vm = the Ising symmetry op¬ 

erators of chain m. That this Hamiltonian is equivalent 
to the Wcluster model with J = 1 and h = 0 can be 
shown by using the following duality transformation 



where now n = 1,..., — 1 and we have defined the 

string operators = Wi^jTf^i+rn- Note that these 
are defined in terms of r Pauli operators! Inverting 
these duality transformations gives the transformations 
given in the main text for the 3-cluster model. 

On the other hand, the A-cluster model with h = 1 
and J = 0 is obtained from by directly identi¬ 
fying and introducing the dual 

operators 


A system of A locally decoupled TFI chains can be 
realised on a single ID system by considering a single 
TFI chain with Ath nearest neighbour interactions only. 
When they are coupled together non-locally such that the 
boundary condition (periodic or anti-periodic) of chain n 
depends on the product ^be Ising symme¬ 

try sectors of all other chains, it has been argued in Ref 
m that the criticality, when all the TFI chains are si¬ 
multaneously critical, changes from the Ising^^ to the 




where now r<^ = j = 0,..., L/(A - 

1) — 1 and n = 0,..., A — 2. 
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Example: 3-cluster model with Ising perturbation 


The spectra of so(A'^)i critical points 


To illustrate the general duality transformations, for 
the 3-cluster model with Ising perturbations they are ex¬ 
plicitly given by 


2 _ y z y 

—^N]^Nj + l^Nj+2^ 

'^3j+2 =^Nj + l<^Nj+2 


and 


X I T r y y z 

'^^3 ~ I \.L^Ni^Ni+l^Ni+2 I Nj'^Nj+l'^Nj+2 

/ 


\i>j 


'3i + l 


n 


y y - \ c, 

^Ni^Ni+l^Ni+2 


'^3j-l-2 —'^Nj+2 I n ^lvi^Wi+l^lvi-l-2 

\i>j 


-z „V „V 


(19) 


Conformal field theory gives a detailed prediction for 
the spectra of one-dimensional critical systems [3] . When 
the finite-size spectrum is rescaled as described in the 
main text, the CFT predicts the energy levels and their 
degeneracies. This information is encoded in the par¬ 
tition function of the CFT, which is composed of the 
partition functions of the left and right moving pieces, 
one for each primary field. In general, for a primary field 
(f>i, with scaling dimension hi, the left moving part of the 
partition function reads 


OO 

E (22) 

ni=0 

(and similar for the right moving part), where the are 
constants, depending on the primary field, and we view 
qi as a formal variable. The total partition function takes 
the form 


(23) 

i 

where the sum runs over all primary fields in the theory. 
In the sector corresponding to primary field (j)i, one can 
expand in powers of q to obtain 


for j = 0, ...,L/3 — 1. Inserting these into (121 with 
N = 3 we obtain the translationally invariant critical 
Hamiltonian 


(1 + Cl,o9i + Co,l9r + C2fiqf + Ci^iqiqr + Co,29r +■■■)■ 

(24) 


Halter = E ^r+Wr+2 + «+l (20) 

i=0 

with a hidden global symmetry described by the 
symmetry operators 

L/ 3-1 

■^0 = n <^Ni<^Ni+l<^m+2^ 

i^O 

L/3-1 

^1= n (21) 

i^O 

L/ 3-1 

1^2 = O'^j + iCr^j+2 

i=0 

that satisfy VqViV 2 = IliE ^i- 


We assumed that the left and right scaling dimensions 
are equal: hi i = h^ i = hi, as is the case for the CFTs we 
consider. The coefficients c„, and the powers 

n of the q variables encode the spectrum as follows: Each 
term in this expansion correspond to degenerate 

states of energy E = 2hi + n. 

The primary field content of so{N)i CFTs depends on 
the parity of N. For odd N the so(iV)i CFTs contain 
three primary fields 1, tp and a with scaling dimensions 
hi = 0, h^ = 1/2 and h^ = N/16, respectively. On the 
other hand, for even N there are four primary fields 1, V'j 
A and A with scaling dimensions hi = 0, /ly, = 1/2 and 
h\ = hx = N/ld), respectively. A thorough derivation of 
the partition functions of so{N)i CFTs can be found in 
Ref. [T5]. Here we just present the predictions for the 
few lowest lying energy levels and their degeneracies that 
are summarized in Tables |T] and [Hj For comparison, we 
also present in Table Ell the CFT prediction for the low 
lying energy levels of the Ising^'^ CFTs. 




FIG. 4. Phase diagrams for general ^-cluster models with 
colour encoding the magnitude of the spectral gap. The cases 
A = 1 and A — 2 correspond to the transverse field Ising and 
XY chains, respectively. For A > 2 The cluster phase occurs 
always around J,h 1, while the (anti-)ferromagnetic and 
the polarized phases occur again for J ^ h,l and h ^ J,l, 
respectively. The fine-tuned incommensurate phases arising 
from the competition between all three terms are gapped mag¬ 
netic phases. 





























9 


N=1 N=3 N=5 



E 

^ni ,n^ 


4*i 

E 

^ni ,n.p 


(pi 

E 

Cn; ,n.p 


0 

1 



0 

1 



0 

1 

1 

1 

0 


1 

1 

6 


1 

1 

20 


2 

2 



2 

27 



2 

160 


1 

1 



1 

9 



1 

25 

Ip 

2 

2 


Ip 

2 

24 


Ip 

2 

150 


3 

3 



3 

88 



3 

785 


1/8 

1 



3/8 

4 



5/8 

16 

(7 

9/8 

2 


<7 

11/8 

24 


a 

13/8 

160 


17/8 

3 



19/8 

84 



21/8 

880 


TABLE I. The lowest lying energy levels E = 2hi+n and their 
degeneracies for the primary field sectors (f>i £ {1, ip, cr} 

of so(A^)i CFTs with N odd. The relevant scaling dimensions 
are hi — 0, h^, = 1/2 and ha = N/16. Degeneracy of Cn,,nr = 
0 denotes no state at this energy. 



(pi 

E 

Cnj ,nr 


0 

1 

1 

1 

2 


2 

9 


1 

4 

Ip 

2 

8 


3 

20 


1/4 

1 

A, A 

5/4 

4 


9/4 

10 


TABLE 11. The lowest lying energy levels E = 2hi + n and 
their degeneracies for the sectors (pi £ {l,i/),A,A} of 

so(A')i CFTs with N even. The relevant scaling dimensions 
are hi = 0, h^ = 1/2 and h\ = h^ = N/16. The fields A 
and A are spectrally indistinguishable, which means that the 
energy levels for these primary fields are degenerate by an 
additional factor of 2. 
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N=1 N=2 N=3 


E 

deg. 

0 

1 

1 

2 

2 

9 

1/8 

2 

9/8 

6 

17/8 

18 

1/4 

1 

5/4 

4 

9/4 

10 


E 

deg. 

0 

1 

1 

3 

2 

15 

1/8 

3 

9/8 

12 

17/8 

48 

1/4 

3 

5/4 

15 

9/4 

54 

3/8 

1 

11/8 

6 

19/8 

21 


E 

deg. 

0 

1 

1 

1 

2 

4 

1/8 

1 

9/8 

2 

17/8 

3 


N=4 N=5 N=6 


E 

deg. 

0 

1 

1 

5 

2 

30 

1/8 

5 

9/8 

30 

17/8 

165 

1/4 

10 

5/4 

70 

9/4 

370 

3/8 

10 

11/8 

80 

19/8 

420 

1/2 

5 

3/2 

45 

5/2 

to 

o 

5/8 

1 

13/8 

10 

21/8 

55 


E 

deg. 

0 

1 

1 

6 

2 

39 

1/8 

6 

9/8 

42 

17/8 

258 

1/4 

15 

5/4 

210 

9/4 

720 

3/8 

20 

11/8 

180 

19/8 

1080 

1/2 

15 

3/2 

150 

5/2 

915 

5/8 

6 

13/8 

66 

21/8 

414 

3/4 

1 

7/4 

12 

11/4 

78 


E 

deg. 

0 

1 

1 

4 

2 

22 

1/8 

4 

9/8 

20 

17/8 

96 

1/4 

6 

5/4 

36 

9/4 

162 

3/8 

4 

11/8 

28 

19/8 

124 

1/2 

1 

3/2 

8 

5/2 

36 


TABLE III. The lowest lying energy levels E < 3 and their degeneracies of the Ising^^ CFTs, for N = 1,..., 6. For each 
Ising^^ CFT the spectrum splits into 3^ primary held sectors, which, for the sake of clarity, are not presented. 



